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ABSTRACT 

Approximately 30% of luminous red giants exhibit a Long Secondary Period (LSP) of variation in 
their light curves, in addition to a shorter primary period of oscillation. The cause of the LSP has 
so far defied explanation: leading possibilities are binarity and a nonradial mode of oscillation. Here, 
large samples of red giants in the Large Magellanic Cloud both with and without LSPs are examined 
for evidence of an 8 or 24 /Ltm mid-IR excess caused by circumstellar dust. It is found that stars with 
LSPs show a significant mid-IR excess compared to stars without LSPs. Furthermore, the near-IR 
J-K color seems unaffected by the presence of the 24 /j,m excess. These findings indicate that LSPs 
cause mass ejection from red giants and that the lost mass and circumstellar dust is most likely in 
either a clumpy or a disk-like configuration. The underlying cause of the LSP and the mass ejection 
remains unknown. 

Subject headings: stars: variables: other — stars: AGB and post-AGB — binaries: general — stars: 
low-mass, brown dwarfs — circumstellar matter — stars: mass loss 



1. INTRODUCTION 

When red giants become more luminous than L ~ 
1000 Lq, they begin to vary with periods of vari- 
ation which fall on six physically distinct period- 
luminosity sequ ences 1 (IWood et al.lll999t llta et afll2004t 
ISoszvhski et all 120071 : Eraser et all I2008D . As well as 
these six sequences for the most luminous red giants, 
there is another sequence (known as sequence-E) extend- 
ing to lower luminosities and which is known to consist 
of clo se binary systems exhibiting ellipsoidal l ight varia- 
tions (|Wood et aLlfT^ISo^ivhski et alJl200l . 

Among the six luminous sequences all but the longest 
in period can be explained by ra dial pulsation in th e fun- 
damental and overtone modes (|Wood et al.l ll999). The 
longest period sequence, sequence-D, consists of stars 
that exhibit variations with a short period, usually cor- 
responding to radial pulsation in a low overtone mode, 
as well as a Long Secondary Period (LSP). It is the LSPs 
that make up sequence-D in the period-luminosity di- 
agram for variable red giants. App roximately 30% of 
all luminous red giants have LSPs ( Wood et all 119991; 
Percy fc Bakos 2001 ISoszvhski et al.l 120071 : iFraser et"aLl 



2008T ) so the LSP phenomenon is very common in the 



late evolution of low mass stars. ISoszvhski et al.l (|2007f ) 
suggest that LSPs can be seen at very low amplitude in 
up to 50% of luminous red giants. The LSPs also s eem 
to occur in red supergiants (jStothers fc Leune]|1971f h 

Since the period of the LSP is approximately 4 times 
longer than the period of the normal-mode radial funda- 
mental mode, the LSP can not be due to normal-mode 
radial pulsation. The most favoured explanations for 
the origin of LSPs are binarity and nonradial g-modes 
(|Wood et al. l fl99llHinkle et al.ll2002t IWood et alF 



Electronic address: wood@mso.anu.edu.au and nicholls@mso.anu.edu. 

1 On these sequences, stars on the first ascent giant branch 
(FGB) have slightly longer periods than those on the asymp- 
totic giant branch (A GB) at the same luminosity llta et al,H2004l ; 
Soszyhski et al. 2007). Some authors designate the corresponding 
FGB and AGB sequences as distinct sequences. 



IDerekas et"aT1l200l ISoszvhskil 120071 IHinkle et al.1l2009T) . 

but these and other e xplanations for the LSPs all have 
signifi cant problems (jWood et al.1 I2004t INicholls et al.1 
2009). For example, stars exhibiting LSPs have veloc- 
ity curves which generally are of similar shape and, in a 
binary model, this would indicate a favoured value for the 
angle of periastron of the orbit rather tha n the expected 
uniform distribution (|Nicholls et al.l [2009T ) . In addition, 
the full amplitudes of the veloci ty curves are close ly con- 
centrated around 3.5 km s _1 ( |Nicholls ~~et al. 2009) which 
implies the unlikely situation wherein all these stars have 
companions of a similar mass of ~0.09 M Q . The problem 
with the nonradial g-mode explanation for the LSPs is 
that g-modes, which have the correct periods, only have 
substantial amplitudes in radiative regions and red gi- 
ants have con vective envelopes w ith only a thin radiative 
layer on top (|Wood et alj 120041 ). Currently, no satisfac- 
tory explanation for the LSPs exists. 

Several of the suggested explanations for the LSPs in 
red giants involve circumstellar dust. The change in vis- 
ible light associated with the LSP is typically of a rather 
irregular nature, and the light variation can be large, up 
to a factor of two. In the binary model, it has been sug- 
gested that the light change could be due to a large dust 
cloud orbiting with t he companion and obscuring the red 
giant once per orbit (|Wood et al.lll999f). Anoth er possi- 
ble explanation for the LSPs (jWood et al .1119991 ) involves 
semi-periodic dust ejection events su ch as those predicted 
by theoretical mod els of AGB stars ([Winters et al.lli"994l 
IHofner et aT1[T995l) . 

Confirmation that dust was involved with the LSP phe- 
nomenon would be obtained if these stars showed a mid- 
infrared flux excess due to absorption of stellar light by 
circu mstellar dust follow ed by re-radiation in the mid -IR. 
a nBoth lHm¥le et al.l ((20021) and lOlivier fc Wood! ((20031) ex- 
amined the mid-IR colors of small samples of stars with 
LSPs in the solar vicinity using IRAS data. They found 
no evidence for a mid-IR excess that would indicate the 
presence of unusual amounts of circumstellar dust. With 
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Fig. 1. — Typical light curves of sequence-D stars ( top 3 panels) 
and o f high mass loss rate stars removed from the IFraser et al.l 
(2008) list of sequence-D stars (lower 3 panels). 

the completion o f the Spitzer Space Telescope SAGE sur - 
vey in the LMC (|Meixner et al.ll2006t iBlum et alj|2006[) , 
it is now possible to search for the presence of a mid-IR 
excess in large samples of stars with LSPs, good con- 
temporary light curves and a known distance. We now 
describe such a search. 

2. SELECTION OF SAMPLES OF VARIABLES 

Th e sample of variab le stars examined here is taken 
from lFraser et al.l (|2008l ) who classified the luminous red 
giant variables in the LMC according to whether they be- 
longed to sequence-D or sequences 1-4 2 . The sequence-D 
stars exhibit light variations at both a primary period, 
which is usually on scqucncc-2 or scqucnce-3, and a Long 
Secondary Period, which falls on sequence-D. Examples 
of the light curves of typical sequence-D stars are shown 
in the top three panels of Fi gure [J 

The sequence-D sample of IFraser et al.l (|2008l) was re- 
fined in a number of ways. To be retained in our sample 
of sequence-D stars, a star was required to show both 
the quasi-periodic LSP as well as a shorter primary pe- 
riod. In an initial examination o f the light cu r ves o f the 
sequence-D stars in the catalog of lFraser et al. ( 2008|) , we 
found stars which have light curves that show a smooth, 
large amplitude variation as well as a variation of mean 
magnitude over very long timcscalcs. Examples of their 
light curves are shown in the lower three panels of Fig- 
ure [TJ These are stars with large-amplitude, Mira-like 
pulsation, high mass loss rates and thick circumstellar 
dust shells and they belong on sequence-1. The dust 
shells cause them to appear fainter in the visible and 
near-IR part of the spectrum than they would appear 
without the dust shell. In the (K,\ogP) plane, they 
thus lie below sequence-1, a nd they can t ake a position 
near or on sequence- D (e.g. lWoodl[2003l ). In order to 
decrease the number of light curves to be examined to 
test for genuine sequence-D characteristics, we restricted 

2 Other authors use sequen ces C, C , B and A in place of se- 
quences 1-4, respecti vely, of IFraser et al.l 1)2008). In addit ion, 
ISoszvhski et al] 1120071 ) find a sh orter period sequen ce that would 
be sequence 5 in the notation of Frascr et al. (2008). 



our sample of sequence-D stars to have a value for the 
LSP falling on the main body of the sequence, defined 
by 2.1 < logP - (14 - if)/4.2 < 2.4. This eliminates 
considerable numbers of the dusty, Mira-li ke stars as- 
signed to sequence-D by IFraser et al.l (|2008[ ) , as well as 
some sequence-E stars. The light curves of all remain- 
ing stars within the main body of sequence-D and with 
K < 12 were inspected and the dusty, Mira-like stars 
and sequence-E stars were removed. We limited our de- 
tailed examination of light curves to K < 12 as it is these 
brighter stars that are used later in the paper for the com- 
parison of mid-IR properties. A small number of other 
stars that did not have the required sequence-D light 
curve shape were also removed: these were normal semi- 
regular variables without obvious LSPs that are typical 
of sequences 2 and 3, and several R Coronae Borealis 
stars. In total, we removed 105 stars from the original 
sequence-D sample. Finally, abou t 5% o f stars appeared 
twice in the lists of IFraser et al.l (|2008f ) under different 
MACHO names (because they were in overlapping parts 
of different MACHO fields): only one entry was retained 
for these stars. 

Our aim is to see if unusual amounts of circumstel- 
lar dust are associated with stars exhibiting LSPs. In 
order to make such a test, we need a sample of stars 
that are identical to the stars with LSPs except that 
they do not have LSPs. This comparable sample is the 
combined group of stars on sequence-2 and sequence-3 
i.e. stars with an oscillation similar to the primary oscil- 
lation period of the LSPVs. Hereinafter, we will refer to 
the variable stars on sequence-D as LSPVs and the stars 
on sequence-2 and sequence-3 as SRVs (semi-regular vari- 
ables). 

Looking at the light curves of the LSPVs and SRVs, 
it is clear that some variations are of low amplitude and 
therefore the periods, an d hence the variabil ity sequence, 
are not well determined. IFraser et al] (|2008f) determined 
periods by two methods - a Fourier analysis and a Super- 
Smoother method. In order to get samples of LSPVs and 
SRVs with reliable sequence determinations, we chose a 
subsample of the LSPVs and SRVs in which the Super- 
Smoother and Fourier period agreed to 5%. We call these 
the P5 samples. The stars in the P5 samples tend to have 
larger amplitudes and more regular variations than the 
remaining stars in the complete samples. The amplitude 
distribution of the two samples is shown later, in Sec- 
tion [3?3l For the stars with K < 12 that are examined 
in detail in this paper, there are 4807 LSPVs and 4568 
SRVs in the complete samples and 1568 LSPVs and 971 
SRVs in the P5 samples. 

3. TESTING FOR CIRCUMSTELLAR DUST USING MID-IR 
FLUXES FROM THE SAGE SURVEY 

Since cool circumstellar dust emits preferentially in 
the mid-IR (longward of ~5 /im) while the stellar photo- 
sphere emits mainly in the near-IR (shortward of ^3 /im) , 
a large ratio of mid-IR flux to near-IR flux can be used 
as an indicator of circumstellar dust. The if- [24] color is 
a sensitive indicator of cool circumstellar dust, especially 
very cool dust with a temperature of ~200K, while the 
if- [8] color would be more sensitive to slightly warmer 
dust with T ~ 400 K. A problem with the if- [24] color is 
that the SAGE survey only detected a 24 fim flux from 
a small fraction of the bright red giants in the LMC. 
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Fig. 2. — if and J-if plotted against If- [8] for the P5 samples 
of LSPVs (small red points, panels 1 and 3 from the bottom) and 
SRVs (small black points, panels 2 and 4 from the bottom). Stars 
on sequence-1 are shown as small grey points. The large amplitude 
Mira variables belonging to sequence-1 are shown as larger points 
that are a darker shade of grey. The magenta crosses are the 10 RV 
Tauri stars in the LMC studied bv lAlcock et al.l (119981 ') that have 
J, if and [8] magnitudes wh ile the blue plus sign s are the 11 RCB 
stars in the LMC studied bv lAlcock et al.l J2001). The lines in the 
bottom two panels are model AGB stars with spherical dust shells 
of increasing optical depth (see text). The top panel shows the 
fraction of the O-rich (J-if < 1.4) SRVs and LSPVs stars brighter 
than if = 12 which lie in if- [8] color bins of width 0.2 magnitudes: 
one-sigma error bars are given. Also shown in the top panel is the 
distribution of sequence-1 stars, both C and O-rich. 

On the other hand, and 8 /im flux was detected for most 
bright red giants. We now discuss both of these options. 

3.1. The K-[8] color 

In Figure [5] we show if and J-if plotted against K-[8] 
for the LSPVs and the SRVs in the P5 samples. Also 
shown are stars on sequence-1: these are fundamental 
mode radial pulsators and they include the Mira vari- 
ables and large amplitude, high mass loss rate, dust en- 
shrouded long period variables. 

It is well known that, in the LMC, red giant stars with 
J-if > 1.4 are mostly carbon star s and those with J-K < 
1.4 are oxygen-rich M or K stars (jCioni fc Habinel l2003h 
The separation of the C- and O-rich stars is seen in the 
(if-[8] , J-if ) diagrams in the lower two panels of Figure^ 
the separation is seen most easily for the grey points of 
sequence-1. These diagrams also show that there is a long 
tail of sequence-1 stars extending from the clump of C 
stars around (if-[8], J-if ) = (1.1,1.8) to redder if- [8] and 
redder J-K colors. This tail contains many C-rich Mira 
variables, defined as sequence-1 variables with a MACHO 
blue amplitude greater than 2.5 magnitudes. However, 
the Mira variables are not confined to this region of the 
(K-[8),J-K) diagram. There are even more of them lower 
in the diagram with J-K < 1.4 and K-[8] < 1.5: these 



are O-rich Mira variables. 

Mira variables are AGB stars of relatively high mass 
loss rate with r oughly spherica l shells, perhaps with some 
dumpiness (see Olofsson 2004 for a review of the evidence 
for dumpiness in Mira winds). In order to get some idea 
of the positions of stars with spherical mass loss shells in 
the (K-[8],J-K) diagram, mo dels of such stars h ave been 
made with the code DUSTY (jlvezic et al.lll999f ) which is 
designed to model the spherical dust envelopes of mass- 
losing AGB stars. Two series of models were made, one 
for C stars and one for O-rich stars. The C star models 
had a central star which was assumed to be a blackbody 
of temperature 2200 K and the grains in the stellar wind 
emanating from the central star were assumed to be 50% 
amorphous carbon and 50% SiC. The O-rich models had 
a central blackbody with T = 2700 K and warm silicate 
dust in the stellar wind. In both sets of models, the 
inner temperature of the dust shell was assumed to be 
1000 K, and the model series had a range of V-band 
optical depths from 0.1 to 10, corresponding to mass loss 
rates of approximately 2xl0~ 7 to 7xl0 -6 M yr _1 for a 
star with L = 5000 L . These are simple approximations 
that do not allow for variation of the stellar temperature, 
inner boundary temperature or grain composition as the 
stars evolve. They are meant to provide a general idea 
of how mass-losing AGB stars evolve in the (if- [8], J-K) 
diagram rather than a precise prediction. 

The upper curve in the lower two panels of Figure [3 
represents the sequence of model C stars with increas- 
ing mass loss rate. This corresponds, within the simple 
model approximations used here, to the tail of sequence- 
1 stars in the upper parts of the (K - [8], J- K ) diagrams. 
This tail thus represents the evolutionary sequence of C 
stars with sp herical mass loss shells o f increasing mass 
loss rate fsee IGroenewegen et al.l [20071 for detailed mass 
loss rate estimates for these stars). The lower curve in 
the bottom two panels of Figure [2] is the equivalent curve 
for O-rich stars. It can be seen that many O-rich Miras 
fall on this curve. 

The position of the high mass loss rate AGB stars 
with spherical shells contrasts with the position of the 
RV Tauri stars shown in Figure [U The latter can have 
very large if- [8] colors but their J- if colors are similar to 
the colors of stars with unobscured photospheres. Since 
RV Tauri stars are luminous post- AGB stars containing 
a dusty circumbinary d isk which radiates in the mid-IR 
(|de Ruvter et al.l [20 051. a low J- if and a large if- [8] is 
likely to indicate a star with a dusty disk. This picture is 
consistent with the J- if colors of the R Coronae Borealis 
(RCB) stars, which are all C-rich. The dust associated 
with them is in neither a disk nor a spherical shell but 
arises from the ejection of puffs of dust in random direc- 
tions, leading t o a highly non-spherical dust distribution 
(|Clavton| [l996K The J-if colors of the RCB stars scat- 
ter between the J-if colors of the RV Tauri stars and the 
high mass loss rate C-rich spherical shell sources, as ex- 
pected in this model since the line of sight to some RCB 
stars would contain dust clouds while in other stars the 
photosphere would be seen unobscured. 

The considerations above allow us to interpret the po- 
sitions of the SRVs and LSPVs in Figure [H The great 
majority of the oxygen-rich SRVs (J-if < 1.4) have if- 
[8] < 0.7, indicating very little mid-IR excess and low 
mass loss rates. Similarly, the C-rich SRVs lie mostly 
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in the domain of the C stars with low mass loss rates 
on sequence- 1. This contrasts with the LSPVs, both C- 
and O-rich, many of which appear to show some evidence 
for a K-[8] excess. For the C stars, the large fraction of 
LSPVs with a #- [8] excess can be clearly seen by compar- 
ing the lower two panels of the figure. The top panel of 
Figure [2] shows the K-[8] color distribution of the O-rich 
SRVs and LSPVs more quantitatively: the larger fraction 
of LSPVs with K-[8] > 0.7 is clearly seen (11% of LSPVs 
have #-[8] > 0.7 compared to 2.3% of SRVs). A two sam- 
ple K-S test gives a probability less than 3.8 x 10~ 7 that 
the SRVs and the LSPVs in the histograms come from 
the same underlying distribution i.e. the higher fraction 
of LSPVs with K-[8] > 0.7 is extremely unlikely to be a 
statistical artifact. Note that in the histograms shown in 
Figure [2J and in subsequent figures, we have restricted 
the samples of stars to those with K < 12 since most of 
the stars that have detectable 8 and 24 (im fluxes have 
K < 12. In addition, first ascent giant branch stars are 
excluded by the restriction K < 12. At 8 fim, we find 
that 99.4% of LSPVs (1561 stars), 99.8% of SRVs (969 
stars) and 99.5% of sequence-1 variables (3964 stars) with 
K < 12 were detected. 

There is one notable feature of the (#-[8],J-#) dia- 
gram for LSPVs. The reddest C-rich LSPVs have J- 
K colors that are similar to those of the unreddened C 
stars, and they are distinctly less than the J-K colors 
of the C-rich circumstellar shell sources in sequence-1 
which have the same K-[8] color. There are also a few 
O-rich LSPVs with large #-[8] colors that fall near the 
RV Tauri stars. This suggests that mass loss from LSPVs 
is not spherically symmetric and that the dust may be in 
a circumstellar disk or some other non-spherical distri- 
bution. These results become much clearer when using 
the 24 /zm measurements from the SAGE survey, as we 
will now show. 

3.2. The #-[24] color 

Figure [3] shows a plot similar to Figure [5] but using 
the if- [24] color. It is immediately obvious from com- 
paring Figures [2] and [3] that only a small fraction of the 
sample was detected at 24 /im. In the P5 samples, for 
objects brighter than K = 12, 27% of LSPVs (424 stars) 
and 28% of SRVs (269 stars) were detected at 24 /im 
whereas essentially all stars were detected at 8 /im. For 
the sequence-1 variables, 52% (2077 stars) brighter than 
# = 12 were detected at 24 /im. 

Comparing the LSPV and SRV samples in the lower 
four panels of Figure El it appears that the LSPVs are 
much more likely to show a #-[24] excess than the SRVs. 
This is confirmed by the histogram of #-[24] colors in 
the top panel of the figure where it can be seen that 
a significantly larger fraction of LSPVs with # < 12 
have #-[24] > 1.5 (72% of LSPVs have #-[24] > 1.5 
compared to 19% of SRVs). A two sample K-S test gives 
a probability less than 2.0 x 10~ 42 that the SRVs and the 
LSPVs in the histograms come from the same underlying 
distribution. 

Given the stronger mid-IR excess exhibited by the 
LSPVs, one would expect a larger fraction of LSPVs to 
be detected at 24 /im than SRVs yet, as noted above, 
similar fractions of LSPVs (27%) and SRVs (28%) were 
detected. This is explained by the inset histogram in the 
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Fig. 3. — The same as Figure l2l but using the if- [24] color rather 
than the if- [8] color for the P5 samples. The dashed line marks the 
24 /im magnitude limit ([24] = 10) of the SAGE survey. The inset 
in the third panel shows the fraction of the P5 sample of LSPVs 
(red histogram) and SRVs (black histogram) detected at 24 fim as 
a function of K magnitude. The histograms of LSPVs and SRVs 
in the top panel include both C and O-rich stars. 

middle panel of Figure [3] which shows the detection frac- 
tion as a function of # magnitude. It is seen that the 
SRVs extend to brighter # magnitudes than the LSPVs 
and that there is almost 100% detection of these brighter 
objects even though they do not have large #-[24] val- 
ues. Most importantly, at all # magnitudes where the 
LSPV and SRV samples both exist, the LSPV detection 
fraction is higher than the SRV detection fraction. 

The C stars and the O-rich stars in the lower two panels 
of Figure [3] fall in distinct groups separated by the line 
J-K — 1.4, the C stars lying above this line. As with 
Figure [2l the J-K colors of the C-rich LSPVs in Figure [3] 
are smaller than those of the C-stars on sequence-1 which 
have high mass loss rates in spherical winds. The picture 
for the O-rich stars is less clear. On the one hand, there is 
a very large group of O-rich LSPVs with J-K < 1.4 but 
red #-[24] colors which fall below the line of DUSTY 
models of increasing mass loss rate in spherical shells. 
This group of O-rich LSPVs seems to follow the sequence 
occupied by the RV Tauri stars with circumstellar disks. 
On the other hand, some sequence-1 Mira variables also 
follow this sequence while others have larger J-K values 
as suggested by the spherically symmetric models. It is 
not clear whether this diversity of behaviour is because 
the winds in Miras are not all spherically symmetric or 
if it is because the models are too simplistic 3 . However, 
the reddest LSPVs with #-[24] > 4 appear similar to the 

3 The O-rich models of Groenewegcn (2006) follow the same 
general behaviour as the DUSTY models shown here but they tend 
to be 0.5 to 1 mag redder in if- [24] for the higher mass loss rates. 
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Fig. 4. — The same as Figure |31 but for stars detected at 24 fim 
in the larger, complete samples. 

RV Tauri stars and they almost certainly have a non- 
spherically symmetric dust distribution. Overall, these 
results indicate that not only do the LSPVs produce more 
dust than equivalent SRVs without LSPs, but this dust 
appears to be in a non-spherical distribution about the 
star, possibly in a circumbinary disk. 

The numbers of stars in the P5 samples are moder- 
ate. However, the results shown above are confirmed by 
the larger number of stars in the complete samples, as 
shown in Figure 3J In the complete samples, for objects 
brighter than if = 12, 31% of LSPVs (1482 stars) and 
28% of SRVs (1285 stars) were detected at 24 /xm. These 
fractions are similar to the detection fractions in the P5 
samples. 

One question that arises is whether the stars that show 
a If- [24] excess are also the stars that show a K-[8] ex- 
cess. We note that there are only 3 stars in the combined 
complete samples with a 24 /im detection but not an 8 fim 
detection, so that it is possible to compare the if- [24] ex- 
cess with the if- [8] excess for essentially all stars. This 
is done in Figure [5] where if- [24] is plotted against if- [8] 
for all stars with if < 12 in the complete samples. There 
are clearly two sequences in Figure [SJ one for C stars and 
one for O-rich stars. Both sequences show a correlated 
increase in if- [8] color with if-[24] color. Thus stars with 
a if- [24] excess do also show a if- [8] excess. Note that 
the LSPVs have a higher fraction of red stars than the 
SRVs (as already demonstrated in Figures [2H1]) • 

3.3. The effect of amplitude on if -[24] color 

It has been shown above that the existence of a Long 
Secondary Period in a red giant seems to induce extra cir- 
cumstellar dust formation. In this situation, it might be 
expected that the infrared excess caused by dust would 
depend on the light amplitude of the LSP. We now ex- 




FlG. 5. — if- [24] plotted against K-[8] for the stars that have 
K < 12 in the complete samples of LSPVs and SRVs. The black 
points are stars with J-K < 1.4 (O-rich star candidates) while the 
grey points are stars with J-K > 1.4 (C star candidates). 
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Fig. 6. — The K-[2A] color plotted against the MACHO blue 
amplitude of the LSP for LSPVs with K < 12 (bottom panel). 
Black dots show stars in the complete sample while stars also in 
the P5 sample are surrounded by red circles. The histogram in the 
top panel shows the fraction of all LSPVs which have K < 12 and 
a detected 24 fim flux in the P5 sample (red line) and the complete 
sample (black dotted line). 



amine this suggestion. 

Figure [5] shows the if- [24] color plotted against the 
MACHO blue amplitude for LSPVs in the P5 and com- 
plete samples. We used the MACHO blue amplitude 
rather than the MACHO red amplitude as the phase cov- 
erage is generally better in blue than in red. 

There is a slight increase of the if- [24] color with 
MACHO blue amplitude. A least squares fit gives 
if - [24] = 0.73Amp(M B ) + 1.47 for the complete sample 
and if - [24] = 0.90 Amp(M B ) + 1-20 for the P5 sample. 
In both cases, the rms scatter of the if- [24] color about 
the fit is 0.50 magnitudes. Given this large scatter, a 
wide range in if- [24] color can be found at any given 
light amplitude. We conclude that the light amplitude 
does not seem to be an important factor in determining 
the amount of mass loss and the mid-IR excess of LSPVs. 

The top panel of Figure [6] indicates that the LSPVs 
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that are in the P5 sample have larger average amplitudes 
than those in the complete sample. This confirms our 
assertion in Section [2] . 

It would be interesting to see if there was any correla- 
tio n between the veloci ty amplitude of the LSP given 
in iNicholls et all (|2009f l and the if- [24] color. Note 
that the radial velocity amplitude associated with the 
LSPs does not show any correlation with light amplitude 
(|Nicholls et a l.l l2009fl. There are 1 6 LSP Vs in the group 
of stars studied bv INicholls elaTl (p009h that have K-[8] 
colors and 6 with if- [24] colors. A plot of the K-[8] and 
if- [24] colors against velocity amplitude for these small 
samples of stars did not show any particular correlation. 

4. SUMMARY AND CONCLUSIONS 

We have shown that luminous red giant stars that ex- 
hibit Long Secondary Periods have larger mid-IR fluxes 
than similar stars without LSPs. This suggests that the 
LSP induces additional mass loss from the red giant, with 
consequent dust formation and an increase in the mid-IR 
flux. The mid-IR flux excess is only weakly dependent 
on the amplitude of the light variation of the LSP. A 
comparison of the near-IR J-if color with the mid-IR 
if- [24] and if- [8] colors indicates that the dust is not in 
a spherically symmetric distribution, and is perhaps in a 
disk. 

Although dust is associated with the LSP phenomenon, 
it is still not clear what causes the LSP in the first place. 
In a binary model, where the velocity curves indica te pre- 
dominantly eccentric orbits (Nicholl s et al.l [2009( 1 . mass 
transfer from the red giant near periastron could lead 
to a circumbinary disk. A non-radial pulsation mode of 
low degree could also produce non-spherical mass loss 
and a non-spherical dust distributi on. The semi-periodic 
dus t ejection events pr edicted by iWinters et al.l (|1994l ) 
and lHofner et al.1 (|1995l ) produce, and indeed require, cir- 
cumstellar dust. However, this phenomenon has not been 
found to occur in theoretical models at the low luminosi- 
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ties where many of the LSPs are observed (even below 
the tip of the FGB). Furthermore, it is also not clear why 
the photospheric velocity should vary with the LSP for 
this purely circumstellar process. 

One additional characteristic of LSPVs is that they 
have a chromosp here which varies with the LSP 
()Wood et al.l 12004b . It is likely that both the chromo- 
sphere and the excess circumstellar dust are manifesta- 
tions of the influence of the LSP phenomenon on mat- 
ter above the stellar photosphere. Magnetic effects are 
a possible source of the chromosphere, but a search for 
magnetic fields in two solar-vicinity LSPVs has put an 
upper limit of 100 Gauss on magnetic fields c overing more 
than about 10% of the surface of these stars (jWood et al.1 
2009) . Currently there is no evidence that magnetic fields 
are the source of the chromosphere. 

In summary, we have shown here that the LSP phe- 
nomenon produces excess circumstellar dust compared to 
stars without LSPs. Unfortunately, since all the postu- 
lated models for the LSP are potentially capable of caus- 
ing mass ejection, the present detection of excess dust 
around LSPVs does not help us distinguish between the 
various models. 
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